Diameter dependence of SiGe nanowire thermal conductivity 
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We theoretically compute the thermal conductivity of SiGe alloy nanowires as a function of 
nanowire diameter, alloy concentration, and temperature, obtaining a satisfactory quantitative 
agreement with experimental results. Our results account for the weaker diameter dependence 
of the thermal conductivity recently observed in Sii-xGccc nanowires (x < 0.1), as compared to 
pure Si nanowires. We also present calculations in the full range of alloy concentrations, < x < 1, 
which may serve as a basis for comparison with future experiments on high alloy concentration 
nanowires. 



The potential interest of nanowires as thermoelectric 
materials has been manifest for more than a decade .^^^ 
Although interest was initially motivated by hopes of tak- 
ing advantage of electron confinement in the structures, it 
soon became clear that another advantage of nanowires 
was their potentially strongly reduced thermal conduc- 
tivity. Advanced techniques enabled the measurement of 
thermal conductivity of single nanowires.^ Large reduc- 
tions in Si nanowire lattice thermal conductivity were 
experimentally reported in 2003, further stimulating re- 
search activities in this area.^ Some astonishingly low 
thermal conductivities have also been recently claimed 
on Si nanowires.'^ 

Many of the initial investigations on thermal conduc- 
tivity reduction by nanostructuring had concentrated on 
ordered crystalline structures. However in recent years, 
the advantages of nanostructured alloys have been under- 
lined for various systems, such as embedded nanodots,!^ 
or nanoporous materials.^ These works show that the in- 
terplay between alloy scattering and scattering by the 
nanostructured features can lead to interesting qualita- 
tive differences between the behavior of the thermal con- 
ductivity n of alloy and non-alloy structures. In particu- 
lar, a slower dependence of n on nano- featured size is ex- 
pected when using alloys. Similar effects may thus take 
place in SiGe nanowires, due to the interplay between 
alloy and boundary scattering. Although the thermal 
conductivity of SiGe nanowires was measured in RefJO^ 
no important reduction was reported there, presumably 
due to the low total Ge concentration. It is only very 
recently that the thermal conductivity of homogeneous 
SiGe nanowires with < x < 0.1 have been reported, 
obtaining remarkable reductions in n below the bulk al- 
loy value. 

In addition to the small reported n values, the exper- 
imental study stressed their much weaker diameter de- 
pendence as compared with pure Si nanowires. Thus, 
the question is whether one can theoretically describe 
and quantify these effects in agreement with the experi- 
mental results. In this Letter, we address this problem, 
and we provide theoretical curves for the dependence of 
n on nanowire diameter, alloy concentration, and tem- 
perature. In addition, we also calculate results for Ge 
concentrations higher than the ones experimentally mea- 
sured so far, allowing for future experimental comparison 



with Ge rich samples. 

A full dispersion theoretical model for the thermal con- 
ductivity of Si nanowires was presented in Ref.l2, Such 
a model only required prior knowledge of bulk Si exper- 
imental thermal conductivity, but did not rely on fitting 
to nanowire measurements. Computing the full phonon 
dispersions may be laborious. However, as it was noted in 
ReffT^^, a very good description of nanowire thermal con- 
ductivity can be obtained without the need to compute 
the full dispersions, just by introducing an additional ad- 
justable parameter, the cutoff frequency cJc, combined 
with a much simpler linear dispersion approximation. 

In this simpler approach, the cutoff frequency is the 
only parameter that relies on nanowire measurements. 
But to fit it, it usually suffices to have measurements of 
i^{T) for just one single nanowire, and one is then able to 
calculate the thermal conductivity of nanowires of other 
diameters without any additional fitting. The simplicity 
of this method makes it very attractive for the analysis of 
other systems. In particular, it can be employed to pre- 
dict the thermal conductivity of SiGe nanowires without 
any prior knowledge of the experimental measurements. 
To do so, only two more things are needed: knowledge 
of the experimental i^{T) for bulk Ge, and knowledge of 
i<i{x) for Sii-a^Ge^^ at a given temperature. 

Basis of this calculation approach are provided in 
Refs. 12l and [HI a summary is given here. The thermal 
conductivity is computed as 

Denoting the Si concentration by x, the different mag- 
nitudes in the above equation are: v^'^ = xv^"^- + (1 — 

^)%Ge^ = ^c,Si{vb/vb,Si), ks = Boltzmauu's con- 
stant, h = reduced Planck's constant, T = temperature, 
and y = hw /ksT . The average speeds of sound for Si and 
Ge are obtained from the experimental sound velocities of 
the transverse and longitudinal acoustic branches, ct and 

^h%{Ge) = (2/3)(cT,5i(Ge))"^ + (l/3)(cL,5i(Ge))"^. 

The cutoff frequency for Si was adjusted in Ref|T2] to be 
about 40 THz. The scattering rate r{uj) is expressed us- 
ing Mathiessen's rule, as a combination of anharmonic 
(r^), alloy (ra), and boundary (rfo) scattering contribu- 
tions: 
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Following the virtual crystal approximation, the an- 
harmonic contribution is approximated as a linear inter- 
polation between Si and Ge: 



(3) 



with r" 



,Si 



BsiUj'^Te^^^^^^ and similarly for Ge. 
Parameters B and C were adjusted to fit the bulk exper- 
imental i^{T) curves (see Table [ll). 



TABLE I: Table of parameters. 

parameter value (unit) 

Vb,s^ 6400 (m/s) 

Vb,Ge 3900 (m/s) 

i^c,s^ 38.8 (THz) 

A 3.01 X 10-^^ (s^) 

Bsz 1.51 X 10-^^ (s/K) 

Bce 2.91 X 10-^^ (s/K) 

Cs^ 139.8 (K) 

Coe 69.34 (K) 
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FIG. 1: (Color online) Thermal conductivity hi vs. tempera- 
ture T for Sii-ccGccc NWs. x is the germanium concentration 
and D is the NW diameter. The lines show calculation results 
using Eq^ the symbols represent experimental results from 
RefHH 

The alloy scattering term is derived as 



Rather good quantitative agreement is obtained for the 
three cases with Ge concentration x = 0, 0.004 and 0.09, 
especially at temperatures above 200K. For the NW with 
X = 0.04 and D — 344nm, the difference between the 
theoretical curve and the experimental data reaches 0.3 
at room temperature. This may be due to the fact that 
K is very sensitive to x when x is small (< 0.05, see Figjs] 
and discussion below), i.e., a small deviation of x can 
make important difference in k. 



100 



0.1 




^ T = 300K 



Q.004 

y x = 0.04 

m x = o.o9 

♦ x = 0(Si) 



10 



100 



1000 



x{l — x)Auj^^ 
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where constant A was adjusted to measurements of 
hz{x) of bulk SiGe alloys for giving the best fitting to 
experimental data.^^ Boundary scattering is included as 
r^^ = V})/D^ where D is the nanowire diameter. The 
model is thus the same as in Ref.8, with the difference 
that the nanoparticle scattering rate in that reference is 
here substituted by the boundary scattering term. Here 
we use the same parameters as in that reference, ex- 
cept for A, which is now fitted by the values of Ref. [131 
rather than those of Abeles.^Thus, no attempt has been 
made to fit the nanowire measurements from Ref. TT via 
adjustable parameters. The very reasonable agreement 
with those results (see below) is quite remarkable, given 
the simplicity of the model and the various approxima- 
tions involved. 

A direct comparison between experimental data and 
our calculated results for four NWs is shown in Figjl] 



FIG. 2: (Color online) Thermal conductivity vs. NW diam- 
eter (log. scale) for Sii-xGex NWs with different germanium 
concentrations at room temperature. The lines show calcu- 
lation results using Eq[l] the symbols represent experimental 
results from Ref .4 11 

We show the dependence of k on nanowire diameter 
m Fig|2l For smah D, K is proportional to whereas 
this dependence becomes slower as D increases. The de- 
viation from linear dependence occurs at smaller D the 
larger the Ge concentration, up to about x ~ 0.5. Obvi- 
ously, as X increases further, the situation reverts, with 
pure Ge nanowires displaying a large linearity range (not 
shown) similarly to pure Si nanowires. This effect is due 
to the coexistence of alloy and boundary scattering con- 
tributions, and it is totally analogous to the effect pre- 
dicted on nanoporous materials in Ref. |9| where the role 
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of thickness was played by the distance between pores. 
The slower dependence is related to the very fast fre- 
quency dependence of alloy scattering. Alloy scattering 
blocks high frequency phonons very effectively, but it is 
totally transparent to low frequency phonons. Thus, the 
thermal conductivity of an alloy is dominated by low fre- 
quency phonons with very long mean free paths, whereas 
in non-alloys contains contributions from a larger range 
of frequencies with shorter mean free paths on average. 
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FIG. 3: (Color online) Thermal conductivity vs. germanium 
concentration for Sii-xGex bulk alloy and NWs with different 
diameters. The lines show calculation results using Eq[l] and 
the symbols represent experimental data from Ref j3l andllll 

Introducing a boundary therefore affects the thermal 
conductivity of an alloy already at rather large values of 
whereas for that same D there is little effect on of a 
non- alloy. In the limit of very small however, bound- 
ary scattering dominates over alloy and anharmonic scat- 
tering at all frequencies. In that regime, of Si and SiGe 
become similar, because the effect of 'bowing' introduced 
by alloy scattering disappears. The competition between 
boundary and alloy scattering can be understood in ana- 



lytical terms from Eq. ([T]) . For high enough T one can ap- 
proximate 1+^, so K (X J^"" (^h/D + A'u 
with A' x{l - x)A. 



For D < VbA 



\ boundary scattering dominates 



at all phonon frequencies, and n (x. D. When D is larger 
than this, however, an upper range of frequencies be- 
comes dominated by alloy scattering. The integral can 
be performed analytically, and it asymptotically tends 
to n (X D^l^ iox D » VbA'-'^uj-^. Thus, the onset of 
the slower D dependence of is directly related to the 
Ge concentration: Donset ^ {x{l — x))~^vi)A~^ijO~^ . The 
aformentioned D^^^ never clearly settles, since it is even- 
tually superseded by anharmonic scattering for larger 
leading to the saturation of n towards its bulk value. 

Finally we plot Kj as a function of x in FigjS) It is found 
that increasing the Ge concentration leads to a very fast 
decrease of thermal conductivity when x < 0.005. Keep- 
ing increasing x, the curves of then tend to saturate at a 
minimum value hZrnin around x = 0.4. Crossing throught 
this minimum points, the thermal conductivities increase 
progressively with the increasing Ge concentration. This 
concentration dependence of NWs is very similar to the 
experimentally-observed one of bulk alloys.^ We note 
that dit X = 0.2 is already close to hZruin- Moreover, it 
can beseen that Kmin of bulk material is at least 7 times 
larger than those of of NWs. 

In conclusion, we have presented the theoretical depen- 
dence of SiGe alloy nanowire thermal conductivity as a 
function of diameter, temperature, and Ge fraction. We 
establish the appearance of a slow diameter dependence 
regime beyond a certain onset diameter, which depends 
on the alloy concentration. The results explain the weak 
diameter dependence reported in a recent experiment, 
and are in reasonably good quantitative agreement with 
those experimental results. Results for higher Ge concen- 
trations beyond the experimentally reported range have 
also been provided, and may allow for further testing of 
the theory by future experiments. 
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